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We show that a quantum system with nonlocal interaction can have bound 
states of unusual type — Isolated States (IS). IS is a bound state that is not in 
correspondence with the S-matrix pole. IS can have a positive as well as a negative 
energy and can be treated as a generalization of the bound states embedded in 
continuum on the case of discrete spectrum states. The formation of IS in the 
spectrum of a quantum system is studied using a simple rank-2 separable potential 
with harmonic oscillator form factors. Some physical applications are discussed, in 
particular, we propose a separable A'^A'^ potential supporting IS that describes the 
deuteron binding energy and the s-wave triplet and singlet scattering phase shifts. 
We use this potential to examine the so-called problem of the three-body bound 
state collapse discussed in literature. We show that the variation of the two-body 
IS energy causes drastic changes of the binding energy and of the spectrum of 
excited states of the three-nucleon system. 



I. INTRODUCTION 



Quantum mechanical bound states are known to have the wave functions de- 
creasing rapidly at large distances r. Usually the bound states are possible at 
negative energies [E < 0) only while at positive energies {E > 0) the system has 
only continuum spectrum states with wave functions oscillating at large distances 
r. Nevertheless von Neumann and Wigner showed long ago 1] that a quantum 
system can have a bound state at positive energy E > 0. Such states are con- 
ventionally refered to as ^continuum bound states^ or as ''bound states embedded 
in continuum^ (BSEC). Von Neumann and Wigner used a local potential in their 
study of BSEC. BSECs are also natural when the interaction is nonlocal (seeQ 
and references therein) or in the case of multichannel scattering (see, e. g., 0,0, Q) . 

A phenomenological nonlocal nucleon-nucleon (NN) potential supporting BSEC 
was suggested by Tabakin 0. The potential predicts the NN data fairly well. 
However it was found that Tabakin potential generates an extremely large binding 
energy for the three-nucleon system 0, B. Such a 'bound state collapse^ was 
investigated by several groups of workers [3, 0, fl^ llsL IT^ with Tabakin and 
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similar nonlocal potentials. All these groups associated the bound state collapse 
with the two-body BSEC and suggested various interpretations of the origin of 
such a puzzling phenomenon. 

In some of these studies |0 0, BSEC was interpreted as an S'-matrix 
pole on the real positive energy axis. This is obviously a mistake since it is well 
known from textbooks |l5j | that the unitarity condition for scattering requires 
\S{E)\ = 1 for real i5 > in the case of elastic scattering {\S{E)\ < 1 for real 
_E > in the case of inelastic scattering), hence the S'-matrix S{E) cannot have 
poles on the real positive energy axis. BSEC appears to be a very interesting state: 
contrary to the conventional discrete spectrum states it is a bound state which is 
not associated with any of the 5-matrix poles p^ . We introduce Isolated States 
(IS) that are, by definition, bound states that do not correspond to the S'-matrix 
poles. Isolated states are a generalization of BSECs: any BSEC is IS, however the 
energy Ej of IS can be also negative, in particular, the ground state of the system 
can be an Isolated State. Generally, if the S-matrix S{E) is known than we obtain 
the energies of the discrete spectrum states by associating the S'-matrix poles at 
negative energies with the discrete spectrum states. The information about the 
IS energy cannot be extracted from the S-matrix, such bound states are isolated 
from the continuum spectrum states. 

In the next Section we present a nonlocal interaction supporting IS. Within the 
J-matrix formalism, it is easy to formulate a realization of this interaction which 
makes it possible to find a simple analytical expression for the S-matrix. This 
interaction is used to study the IS formation when the Hamiltonian parameters 
are varied, the IS contribution to the Levinson theorem, etc. 

We address the bound state collapse problem in Section III. We fit the param- 
eters of our exactly solvable nonlocal interaction model supporting IS to the NN 
scattering data. The IS energy Ej is arbitrary because it is not related to the 
S'-matrix. We show that varying the IS energy Ej (note that the S'-matrix is 
unaffected by this variation) we produce great changes of the three-body binding 
energy: some Ej values bring us to the system with three bound states with 
extremely large (few GeV) binding energies; with larger Ej values we obtain the 
•^H nucleus with two bound states (the binding energy of the ground state is of 
the order of few hundreds MeV); the further increase of Ej results in the further 
decrease with Ej of the binding energy of the '^H nucleus which has a single bound 
state; if the IS energy Ej is large enough the '^H nucleus becomes unbound. The 
experimental value of the "^H binding energy can be exactly reproduced if some 
particular positive energy Ej is taken when IS appears to be BSEC. Therefore 
the problem of the three-body bound state collapse does not exists as a general 
problem: this problem arose only due to the use of very restricted models of inter- 
action supporting BSEC for the construction of NN potentials, i. e. this problem 
is inherent for such interaction models only. 

We present in Section IV a short discussion and compare our results with the 
results of other authors who studied the three-body bound state collapse problem. 
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II. SIMPLE POTENTIAL MODEL SUPPORTING ISOLATED STATE 

The radial wave function ^''£;(r') satisfies the Schrodinger equation 

{T^ + - E)¥Eir) ^0, (1) 

where E is the energy, I is the angular momentum, T' is the kinetic energy operator 
and is the potential energy. 

Let i = 0, f, 2, ... be a complete basis. The Hamiltonian matrix 

Hij = is generally infinite and the wave function ^'^(r) at any energy E 

can be expressed generally as an infinite expansion in basis functions, 

oo 

*kM = E"^(^)N)- (2) 

4=0 

However, at some particular energy Ei the infinite Hamiltonian matrix Hij can 
have a finite eigenvector, and the wave function \l//(r) at this energy is expressed 
as a finite expansion in basis functions, 

M 

*z(r)=E (3) 

i=0 

The wave function ^'/(r) rapidly decreases with distance r since it is a superposi- 
tion of a finite number of functions. Therefore at energy Ej we have a bound 
state of a particular type hereafter refered to as an Isolated State (IS). 

Clearly, we have the IS solutions of the type ^ in the case when the Hamiltonian 
matrix Hij is block-diagonal, 

H.,=H^P®Hg\ (4) 

where the (AI -f- 1) x (M -I- 1) submatrix H^p is defined in a finite-dimensional 
subspace spanned by the basis functions \i) with i = 0, 1, ... , M. The infinitc- 
dimensional submatrix H^j is defined in the orthogonal supplement to this sub- 
space. Any eigenvector of the submatrix Hij gives rise to the wave function of the 
type lj2Jl, i. e. each of the submatrix Hij eigenvectors is associated with IS. 

The scattering state wave functions with oscillating asymptotics at energies 
E > 0, can be expressed only as a superposition of an infinite number of 
functions, 

oc 
i=M+l 

where are eigenvectors of the infinite submatrix H^j\ The iS-matrix 

and scattering phase shifts are defined through the asymptotics of the functions 
©, hence they are governed by the structure of the submatrix H>,' only. The 
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energy Ej of IS and its other features are dictated by the structure of the other 



submatrix H^p that is generally independent from H^^'. Therefore we cannot 
expect that the S'-matrix has a pole at the IS energy Ej. We can define the 
Isolated State as the bound state that is not associated with any of the S'-matrix 
poles. The IS energy Ei can be positive, and in this case it appears to be the 
so-called bound state embedded in continuum (BSEC). Any BSEC is IS. The IS 
energy Ej can be also negative, in particular, the ground state of the system can 
be isolated. Thus IS can be treated as a generalization of BSEC on the case of 
arbitrary (negative or positive) energy. 

The J-matrix formalism 17] makes it possible to study IS properties and to 
formulate a simple exactly-solvable model of a system possessing IS. In this con- 
tribution we use the oscillator basis; the exactly-solvable model of IS can be also 
easily formulated by means of the J-matrix formalism with the Laguerre basis. 

The idea of the block-diagonal structure of the Hamiltonian matrix Q can be 
easily realized if the interaction between the particles is described by a separable 
nonlocal potential of the rank N + 1, 



r(2) 



N 



(6) 



with the harmonic oscillator form factors 

^ni{r) = {-IT 



ror(n + Z + |) 



i+i 



exp 



'2^ ) ^" \Z2 



' 



(7) 



Here = (h/mcu) is the oscillator radius and L"{x) is the Laguerre polynomial. 

In the J-matrix method, the wave function has a form of series in terms of 
functions 0, 



(8) 



The coefficients Xn{E) for n > N are given by the formula 
X,i{E) = Sni{p) COS Si + Crd{p) sinSi , 



(9) 



where p ~ ^2E jfiLo is the momentum, Sni (p) and Cni (p) are the eigenvectors of 
the infinite tridiagonal matrix of the kinetic energy r^„/ . The following analytical 
expressions 0| can be used to calculate Sni {p) and Cni (p) '■ 



Snlip) = 
Cnlip) = 



2T{n + l 



3m 2 
2 



r(n-f 1) 
2r(n+ 1) 



P 



P 

p^.;_^ 3^exp(-y) / + ^; p^) , 



(-1)' 



, , -T— exp( n- 



(10) 



-l + -;p^ 



(11) 
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H(2) 




Fig. 1: The structure of the Hamiltonian matrix. 



The phase shift Si in the partial wave with the angular momentum I can be cal- 
culated as 

, r _ Sni{p) - PNn{E) Sn+iAp) /-,r,^ 

' Cm{p)-PNN{E)CN+Upy ^ ' 



where 



Pn 



(^) = -Er^^''-'+i' (13) 



(n = 0, 1, ■■■,N) are the eigenvectors and are the corresponding eigenvalues 

of the truncated Hamiltonian matrix 11^^' — ^/m' + Kin' i^' = 0, 1, A''). The 
coefficients Xn{E) ior n < N can be found by the formula 

Xr,{E) = PnN{E)XN+l{E). (14) 

We are considering the case when the Hamiltonian matrix is block-diagonal. We 
note that the kinetic energy matrix in the oscillator basis is tridiagonal. Hence with 
the interaction © we can obtain the Hamiltonian matrix in the oscillator basis of 
the type I0J that has the structure shown in Fig.^ Solid lines schematically show 
the infinite tridiagonal kinetic energy tail of the Hamiltonian matrix, the rest non- 
zero matrix elements arc displayed by two boxes representing submatrices H^^^ 
and _ff (we include the infinite tridiagonal kinetic energy tail in the submatrix 

Obviously, the eigenvectors of the submatrix iJ^^^ are the eigenvectors of the en- 
tire Hamiltonian matrix H, too. The corresponding wave functions are decreasing 
with r similarly to the bound state wave functions since they are superpositions 
(PJ of a finite number of the oscillator functions \i) — </?„;(r) given by Eq. {Tj). 
However these wave functions have an unusual asymptotics exp[— r'^/(2rp)] in- 
stead of the standard one ~exp(— ar). All eigenstates of the submatrix H^^^ are 
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ISes. It is obvious that generally the submatrix iJ*^^^ may have positive or/and 
negative eigenvalues £^^\ If e^^^ > 0, the corresponding state appears to be BSEC. 
If e'j}^ < 0, the corresponding state appears to be a bound state of a specific type. 

The continuum spectrum states with the oscillating asymptotically wave func- 
tions as well as the conventional bound states with the e~"'"-type wave function 
asymptotics, can be expressed only as infinite series (0) of the oscillator functions 
\i) ~ ifniir) given by Eq. {Tj). They are generated by the submatrix H^^\ The 
scattering and conventional bound state wave functions Q are obviously orthogo- 
nal to the IS wave functions JSJ- Thus the scattering and usual bound state wave 
functions have node(s) at a small distance (~ro) in the presence of IS(es) similarly 
to the wave functions generated by the potentials with forbidden states [l^l • 

The asymptotic behavior of the scattering state wave functions is characterized 
completely by the 5'-matrix Q . Hence the structure of the 5- matrix is governed by 
the infinite-dimensional submatrix iJ^^^ that is independent from the submatrix 
The IS energies s^i}\ on the other hand, are controled by the submatrix 
and are independent from the submatrix H'^^\ Varying matrix elements 
^nn' '^^ submatrix H^^'> one causes variation of the IS energies e^^"* without 
affecting the 5'-matrix. Thus, the energy of IS ej^^'' is not in correspondence with 
the location of the S'-matrix poles. Using symmetry properties of the 5'-matrix as 
a function of the complex momentum p 3], it is easy to show 16] that the energy 
of BSEC is not in correspondence with any of the S'-matrix poles. An interesting 
new point, so far as we know never discussed in literature, is the appearance of 
the discrete spectrum states, i. e., of ISes with negative energy e^^"* < 0, that are 
divorced from the S'-matrix poles. So, IS being the state with the asymptotically 
decreasing wave function and with the energy at which the S-matrix does not have 
a pole, can been treated as a generalization of BSEC on the case of the discrete 
spectrum states. 

We examine in more detail the formation of IS in the spectrum of a quantum 
system with nonlocal interaction using as an example a simple analytically solvable 
model. The simplest realization of the situation depicted in Fig.^ corresponds to 
the case when the submatrix i?^^^ is a 1 x 1 matrix and the separable potential 
© is of the rank 2, i. e. N = 1. In this case, IS arises due to the cancellation 
of the potential energy matrix elements Vbi = Vio and the kinetic energy matrix 
elements Tqi = Tio = — Vbi that results in Hqi = -ffio = 0. The IS energy Eq is 
equal to the diagonal matrix element -ffooi = ^^oo- It should be stressed that 
-ffoo can take an arbitrary value in our model. Using Eqs. H10II13|I we obtain the 
following expression for the phase shift Si : 

Soiip) {[Uii(£o ~E)-f3] (T^o -E) + {T},,)' (eo - E)\ 
tan (5/ = p ^ ■ 

Coiip){[Vnieo - E) ~ f3] i^, - E) + (T^,) (e„ - E)} - "^''^^XiP) 

(15) 

where (3 = H^^ . 

Suppose eo > and Vu > 0. The evolution of the s wave phase shift So{p) as 
P = Hq^ = H^Q tends to zero is shown in Fig. |21 If /9 7^ there is a resonance at 
the energy E Eq oi the width F that decreases as /3 is reduced. When (3 = 0, 



7 




-180^ 



Fig. 2: The evolution of the phase shift 5o(p) as /3 = Hqi — H^o tends to zero. Dashed, 
dotted, and sohd curves are the phase shifts obtained with P = Pi, P2 and Pz, respectively; 



BSEC arises as the resonance of the zero width producing the jump of the height 
TT of the phase shift Sq at the energy E = Eq. This spurious jump should be 
ehminated that results in the 5q{Q) increase by an extra tt if we suppose, as usual, 
that (So (00) = 0. Thus when applying the Levinson theorem Q to the system 
pertaining IS, one should count IS as a usual discrete spectrum bound state. Such 
behavior of the phase shift is typical for systems pertaining BSEC that have been 
studied in various models 0, IigL . Thus our model represents an alternative 
simple analytical approach in the study of BSEC. 

As for the S-matrix, it is given by the following expression: 



C. 



Si 



(-) 

0/ 



(p) { [Vl,{s, -E)-(3] {U, -E) + {so E)} 



p[Vl,ieo-E)-l3] 



Cin\p) { [Vhie, -E)-f3] {T^, - E) + [T^.f {so - E)} 



p[Vl^{eo-E)-l3] ' 

(16) 

where c'^\p) — Cniip) ± iSniip)- The single S'-matrix pole on the unphysical 
sheet tends to the real energy E = Eq a,s (3 tends to zero. However in the limit 
P = 0, the factor {eq — E) in the numerator of Eq. (|15|l cancels the same factor in 
the denominator and the singularity at the energy E = Eq disappears: 



S, = 







-E)A 




■^Soiip) 










■^Soi(p) 



(17) 



This illustrates the mechanism of the 5-matrix pole loss in the limit /3 — > when 
the resonance transforms into BSEC. The nontrivial result is that if IS is a bound 
state (eo < 0), it does not generate the S'-matrix pole, too. 
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The above results can be easily generalized as the following statement. 

Let all eigenvalues of the truncated Hamiltonian matrix he non-degenerate. 
Then Isolated States occur in the spectrum of the system with the rank-{N + 1) 
separable interaction ^Bjl if and only if the truncated matrix and its principal 
minor H^~^ of the rank [N — 1) have common eigenvalues. The number v of the 
common eigenvalues is equal to the number of the Isolated States. These eigen- 
values and the corresponding eigenf unctions are just the energies and the wave 
functions of the Isolated States. 

The equivalent formulation is: 

The system with the nonlocal interaction 0) has the Isolated State at the energy 
if and only if is the eigenvalue of the truncated Hamiltonian matrix and 
the corresponding eigenvector U'^ has the last component = 0. 



III. PHENOMENOLOGICAL iViV POTENTIAL WITH ISOLATED 
STATE AND THE THREE-NUCLEON SYSTEM 

The simplest rank-2 separable potential (jHl supporting IS discussed in the 
previous Section, was used to fit the NN singlet ^sq and triplet "^si scattering 
phase shifts. The oscillator function parameter huj = 500 MeV. To ensure the 
existence of IS, we should set the off-diagonal potential energy matrix elements 
Vol = Vio — —TiQ. The phase shifts are independent of the matrix element Vqq 
that governs the IS energy. The only parameter responsible for the phase shifts is 
the matrix element Vn. 

The singlet phase shifts obtained with Vf^ = —0.7315 huj and the triplet phase 
shifts obtained with Vii — — 0.81512 fiij are shown in Figs. 13 and 0] respectively. 
These values of Vfi and Vii are seen to reproduce with a reasonable accuracy the 
scattering data. 

The deuteron ground state energy should be obtained by the calculation of the 
5'- matrix pole as is discussed in Refs. |23,|^- Since in our case the S'-matrix is 
given by the expression IjlTI) , the ^-matrix poles can be calculated by solving the 
equation 



=0. (18) 



The bound state (deuteron) should be searched for in the triplet "^si wave, i. e. 
I = and should be used as Vfi in Eq. 118|l . We are searching for negative Ed 
value and pd — i\/2\Ed\lhhj fitting Eq. (|18|l . The deuteron wave function can be 
calculated using the J-matrix formalism (see also the discussion in Refs. [23, l2l| 
and references therein). Our very simple potential with the only fitting parameter 
Vli provides a good description of the deuteron energy Ed — —2.22496 MeV 
and rms radius \J (r^) — 1.87 fm (the respective experimental values are E"^^"^ = 

-I.Tl^hlh MeV and ^{fl^^) = 1.9676 fm). 

The obtained singlet and triplet s wave NN potentials are used in the calculation 
of the triton bound states. As in the other studies of the three-body bound state 
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collapse 01111111111111111 El, we do not allow for the interaction in other 
partial waves. We perform a conventional variational calculation of the S = T = 
J = \ three-nucleon states with the three-body oscillator basis allowing for all 
components with the total number of oscillator quanta N = 2n + I < 32. The 
energies of the ground state and of the lowest excited states are shown in Fig.Ofor 
different Ej values (we are varying both the triplet Vqq and the singlet Vqq potential 
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Fig. 5; The triton ground and excited state energies vs the IS energy Ej. 

energy matrix elements so that Vqq = Vqq ; in other words the IS energy in the 
triplet state is equal to the IS energy Ef in the singlet state: Ei = Ef = Ej). It 
is seen that the variation of the IS energy Ej (that does not affect the phase shifts 
and the deuteron properties) results in the drastic changes of the triton binding 
energy and of the spectrum of excited S = T = J = ^ states. When the IS 
energy Ej is small enough, the three-nucleon system collapses (the binding energy 
becomes extremely large) ; two excited states are bound in addition to the ground 
state. The triton ground state energy increases with the IS energy Ej; the same 
is true for the energies of the excited states. At some Ej value the second excited 
state becomes unbound; at some larger Ei value the first excited state becomes 
unbound, too; however the triton binding energy is still too large. Nevertheless 
the further increase of Ej results in the increase of the triton ground state energy 
and at some Ej value the three-nucleon system becomes unbound. 

We see that we really have the three-body bound state collapse for some values of 
the IS energy Ej. However the collapse disappears for larger Ej values. Therefore 
our general conclusion is that the three-body bound state collapse problem does 
not exist as a general problem. In the previous studies of the three-body bound 
state collapse, very restricted potential models were used that did not allow to 
change the BSEC energy i?BSEC- The i?BSEC value used in these studies caused 
the extreme overbinding of the trinucleon. This is clearly the problem of the 
particular potentials used and not the general problem of the three-body bound 
state collapse inherent for all two-body potentials supporting BSEC. In our model 
any trinucleon binding energy can be obtained by phase equivalent variation of the 
IS energy Ej. In particular, we can fit the Ej value to the triton binding energy. 
Setting Ej = Ef = Ej = 189.525 MeV, we obtain the triton binding energy of 
Et = 8.47307 MeV in our 32huj variational calculations. The accuracy of the 
variational approach is improved if the three-body ^-matrix pole is calculated as 
is discussed in Ref. [23|: in this case the triton binding energy of Et — 8.4748 MeV 
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is obtained. 



IV. DISCUSSION 

We suppose that the most important feature of BSEC is that this bound state is 
not in correspondence with the S'-matrix pole. We introduce IS as a generahzation 
of BSEC on the case of arbitrary (positive or negative) energy. We propose a simple 
model of non-local interaction elucidating the formation of IS. To investigate in 
detail the formation and the features of IS, it is natural to employ the J-matrix 
formalism that makes it possible to obtain analytical expressions for the S'-matrix 
and other observables, i. e. to formulate the exactly solvable model of IS. 

The general results given above can be straightforwardly extended on the cou- 
pled channel case, on the case of the separable finite-rank nonlocal potentials with 
the Laguerre form factors , etc. 

It is interesting that IS naturally appeares in the standard nuclear shell model. 
It is well-known (the Dubovoy-Flores theorem) J2, 23] that the lowest eigenvalue 
Eq (w) obtained in the shell model calculation in the Nhio model space, coincides 
with the lowest eigenvalue e^~^(w) obtained in the shell model calculation in the 
(N — 2)hLU model space if the parameter Hu! of the oscillator basis minimizes the 
eigenvalue £^~^(ci;). According to the statement formulated in Section III, the 
ground state obtained in the conventional variational nuclear shell-model calcula- 
tions appears to be IS. 

Using our exactly solvable model of IS, we obtain a simple nonlocal NN potential 
describing the singlet and triplet s wave phase shifts and deuteron properties. 
With this NN interaction we study the triton properties and examine the three- 
body bound state collapse problem. We can vary the IS energy Ej in our model 
without violating the description of the phase shift and deuteron properties, i. e. 
without altering the on-shell properties of the interaction. However the off-shell 
properties of the interaction are modified when Ej is varied and this is clearly seen 
in the strong i?/-dependence of the triton ground state energy Et . The three-body 
system collapses at some Ej values and does not at other Ej values. We can even 
fit the position of IS to the triton binding energy. Therefore we conclude that the 
discussion of the three-body bound state collapse in literature HollTH IT^ IT^ IT^ 
arose only due to the use of unsuccessful NN potential supporting BSEC; other 
two-body interactions supporting BSEC provide a correct description of the three- 
nucleon binding. 

Various interpretations of the nature of the three-body bound state collapse 
have been suggested in literature ^ J^, J^J . We cannot accept most of 

them. 

For examples, Nakaichi-Maeda and Pantis et al 0| supposed that the three- 
body bound state collapse arises due to the nodel behavior of the deuteron wave 
function in the case of potentials supporting BSEC. Delfino et al supposed 
that collapse originates from the structure of the BSEC wave function which is 
identical with that of the Pauli forbidden state. In our model, we shifted the IS 
energy changing neither the deuteron wave function nor the IS wave function. We 
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show that the three-body bound state collapse can be eliminated by the increase 
of the IS energy only, therefore we cannot agree that the collapse originates from 
the structure of the deuteron or IS wave function. 

Delfino et al ^if supposed that the three-body bound state collapse is a man- 
ifestation of the Thomas effect (the increase of the triton binding energy when 
the range of the NN potential tends to zero) . The range of the general nonlocal 
interaction cannot be established unambiguosly. Delfino et al suggested to use 
the average kinetic energy (T) of the two-body bound state as an indicator of the 
potential range and demonstrated that BSEC is accompanied by an increase of 
(T). We cannot agree with this conclusion. In our model the bound state wave 
functions are unaltered when the IS energy is varied, hence the variation of IS 
energy does not cause changes of (T) of any of the bound state. Therefore with 
the same values of (T) of two-body bound states, we obtain either the collapsed 
triton at some Ej values or normally-bound (or just unbound) triton at other Ej 
values. 

Rupp et al |10| note that the triton binding energy increases when the BSEC 
energy grows. We obtain an opposite result in our investigation: the triton binding 
energy decreases with the BSEC energy. We note here that in our studies we shift 
BSEC phase equivalently while Rupp et al base their conlusion on the results 
obtained with a number of potentials providing different phase shifts. In other 
words, they use potentials with different on-shell properties that can mask any 
effect of the BSEC energy variation. 

It is interesting that it is easy to change the BSEC (or any other bound state) 
energy phase equivalently. The BSEC wave function ^'bsec fits the Schrodinger 
equation 

-ff^BSEC = -E^bsec'I'bsec- (19) 
Let us define a new Hamiltonian 

H' = if + A|^'BSEC)(*BSEC|. (20) 

The wave functions '5^(7') of all the rest (bound or continuum spectrum) eigen- 
states of the Hamiltonian H, fit the Schrodinger equation 

H^Eir) ^ E^Eir). (21) 

Clearly ^'^(r) will be also the eigenfunctions of the Hamiltonian H' since 5'bsec 
is orthogonal to any \E''g(r). At the same time, the BSEC energy i?BSEC will be 
increased by A. We are sure that if Rupp et al employed this method of varying 
the BSEC energy, they would obtain the natural result that the triton binding 
energy decreases with E'bsec- 

The Hamiltonian (|2(JII can be defined with the projector I'^sJi'^Eb I '-'^ other 
bound state ^'^^(r) at the energy E^. In this case we obtain another simple model 
of IS. Really, none of the continuum spectrum wave functions ^'^(t') is altered by 
the projector (^'^^^ |- Hence the Hamiltonians H and H' provide the same 

S'-matrix. Thus the projector Aj^''^;^) (^''^j^ | increases the energy Ei, of the bound 
state by A but does not shift the S'-matrix pole at the negative energy Eb. As a 
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result the conventional bound state at the energy Ef, is transformed into IS. We 
note here that the Hamiltonians with the projection operators A||l/pJ(^'^J are 
widely used, e. g., to project out Pauh forbidden states (see, e. g., ^4] ) . 

We do not suppose that it is needed to criticize other interpretations of the origin 
of the three-body bound state collapse suggested in various papers. Our general 
conclusion is that the collapse is inherent for very restricted models of potentials 
supporting BSEC only and does not appear if other two-body interactions with 
BSEC are used. Therefore for the proper interpretation of the three-body bound 
state collapse origin, one should carefully study this restriction and reveal what is 
wrong with it. Up to the best of our knowledge, it was never done. 

This work was supported in part by the State Program "Russian Universities" 
and by the Russian Foundation of Basic research. Grant No 02-02-17316. 
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